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raa 1*4 ep (5) . x C5) . y (5) . z (5) . xs (5) . ys (5) , zs (5) 
rea 1*4 ai (200) . bl (200) . oi (201) 



234556 
fO=0. 05 

dlh=15.0/fO ! cm 
na^16 
write 10) 

formate Enter the location of x. y.z (om) 
read(*. *. end=90) xp. yp, zp 



(3) 



•$) (2) 



call marray (xp. yp.zp.na, cl) 
do i=2, na+1 

verr^rand i j) 

al (i-1)=cl (I) ♦ (1, (K(verr-0. 5)*0. 001)-cl (1) 
end do 

writeC*.*) ' ^L(cm)\ (al (i) J=l,na) 
writeC*,*) 

call mcycleCna, dih, al). (5) 

do 1=1.5 

ep(i)=1.0e20 
end do 

do ix=-30,3O 

xp=float(ix)*10. 0 
do ty=-"30, 30 

yp=float(iy)*lQ.O 

do lz=--30.30 



\ noise 0. 1 % 



] 



(4) 



zp=float(iz)*10. 0 



(6) 



(7) 



call fliarrayCxp. yp, zp. na, ol) 
do i=^2. na+1 

bl(i-!)=cl(i)-cl(l)-al<M) (8) 
and do 

08 1 1 moycle(na.d lh.bl) (9) 
ar=0. 0 
do i=1.na 

er=8r+bl Ci)**2 (10) 
end do 
do i=1.5 

if (er . It. ep(i)) then 
if (I .n«. 5) then 
do j=5. 

ep(j)=6p(j-1) 
x(j)=xa-1) 
y(j)=y(j-1) 
2(j)=z(j-r1) 
end do 
end if 

ap(l)»er 
x(i)=xp 

y U>=yp 
£(i)szp 
go to 30 
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end if 
end do 
cont i nue 
end do 
end do 
end do 

do i=1.5 



end do 

writeK*) ' RMS error (cm) x y 

do i=n5 

write(*,*) sqrt(ep(i)/float(na)).x(i).y(i).z(0 
and do 

do 

xO=x8(m) 
yOsys On) 
20=zsOft) 
do ix=-15J5 

xp=f loat (ix)+xO 
do iy=-15. 15 

yp=f loat(iy)+yO 
do iz=-15, IS 

2p=f loat <i2)+20 



call narrayCxp.yp, zp. na.cl) 
do is2. na4-1 

bl(iH)=cl(i)-^l(1)-al(i-l) 
end do 

cali mcycleCna.dlh. bl) 
er=0. 0 
do 1=1. ra 

er=er+bl {i)**2 
end do 

do 1=1.5 

if (er . It. ep(l)) then 
If (i .ne. 5) then 
do j=5. i+l.-l 
ep (j)=Bp(j'-1) 





end do 



end if 



ep(i)aer 
x(0=xp 




go to 35 



end if 



end do 
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continue 
end do 
end do 
end do 

end do 
write (5^, *) 

write(*,*) sqrt(ep(1)/float(na)).x(1),ya),2(l) (U) 

write(*, *) 
go to 20 

stop 
end 

subroutine marray (xp, yp. zp, na, cl) 
real*4 cl (1) 

o I (1 ) =sqr t (xp5^xp+yp*yp+ (zp^SO. 0) ^2) 
do i=2. na+1 

ixx=i/3 

lyy=i-ixx*3 

xm=f I oat ( i xx-1 ) *50. 0-1 0. 0 
ym=f I oat ( i yy-1 ) ♦50. 0+1 0- 0 
o I < I ) =sqrt ( (xp-xm) 'W«2+ (y p-ym) 'f*2+zp*zp) 
end do 

return 
end 

subroutine mcyoteCna, dih, at) 
real^ a I (1) 

do i^t.na 
continue 

if (al (i) .ft. dih) then 

al(i)=aT(i)"d!h 

if (al(i) . ie. dIh} go to 46 

go to 40 
end if 
continue 

if (aUi) , It. -dIh) then 

al(i)=al(i)-Kdth 

if (aid) .ge. -dIh) go to 46 

go to 45 
end if 
continue 
end do 

return 
end 
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Enter the location of x, y, z (cm) : 152,-203.56 
^L(cm) 67.67562 -38.21133 -1.487458 
-69. 24731 -27. 88023 16. 30007 -91. 74537 

0.9732714 -102.0754 -54.30361 -5.570741 
-51.46763 -3.269386 



39. 09471 

-46. 11990 
-98. 28325 



RMS error 
0. 6834297 
0. 8562734 
1. 116775 
I. 163736 
1.216863 



(cm) 



150.0000 
150.0000 
150.0000 
160.0000 
160.0000 



8.4395386E-02 152.0000 



-200. 0000 
-190. 0000 
-200. 0000 
-230. 0000 
-220. 0000 

-203. 0000 



60. 00000 
50. 00000 
50. 00000 
70. 00000 
60. 00000 

56. 00000 



Enter the location of x, y. z (cm) : 22,123.-89 
A-im) 5.506481 57.46710 16.50204 -17.27929 

55.74849 14.06553 -20 41722 66.89948 28.19106 

-2.332703 89.04320 55.22502 29.83902 119.4193 

90.37129 69.39222 



RMS error 
1.445567 
1.754374 

1. 951296 

2. 345274 
2. 709345 



(cm) 



20. 00000 
20. 00000 
20. 00000 
20. 00000 
20. 00000 



6. 2024966E-02 22. 00000 



130.0000 
130.0000 
120.0000 
120.0000 
140.0000 

123.0000 



-90. OOOOO 
-100. 0000 
-80. OOOOO 
-90. OOOOO 
-100. 0000 

-89. OOOOO 



Enter the location of x. y, z (cm) : 60.161,5 
^Kcm) -23.45399 32.54938 -13.85323 

21.66080 -27.96993 -77.36571 

-74. 96463 36. 05470 -9. 367880 

18.86572 -15.62937 



-57.41031 
22. 85288 -26. 38201 

51.50449 59.00156 



RMS error 
1.358104 
1.400364 
1.561480 
1.779230 
1.850774 



(cm) 



60. OOOOO 
60. OOOOO 
60. OOOOO 
60. OOOOO 
60. OOOOO 



160. 0000 
160. 0000 
170.0000 
170.0000 
150. 0000 



4. 4650473E-02 60. OOOOO 161. 0000 

Enter the location of x. y, z (cm) : 



10. OOOOO 
0. OOOOOOOE+00 
0. OOOOOOOE+00 

10. OOOOO 

10. OOOOO 

5. 000000 
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